Navier-Stokes Equation by Stochastic Variational Method 
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We show for the first time tliat the stochastic variational method can naturally derive the Navier- 
Stokes equation starting from the action of ideal fluid. In the frame work of the stochastic variational 
method, the dynamical variables are extended to stochastic quantities. Then the effect of dissipation 
is realized as the direct consequence of the fluctuation-dissipation theorem. The present result reveals 
the potential availability of this approach to describe more general dissipative processes. 



I. INTRODUCTION 

In almost all branches in physics, the variational prin- 
ciple is one of the important guiding principles, and also 
serves as a powerful tool in practice In particular, this 
principle is indispensable method to deal with the sym- 
metry property of a system, for example, the relativistic 
covariance. Once the symmetry of a system is expressed 
as the invariance of the action under the transformations 
of the symmetry in question, then the resultant equations 
of motion automatically satisfies this symmetry. 

The usual variational principle, however, is not ap- 
plicable when irreversible dynamics is present. Because 
dissipation involves energy exchange processes between 
macroscopic and microscopic motions, associated with 
the entropy production mechanism. These microscopic 
degrees of freedom are not included as the system de- 
grees of freedom in the Lagrangian. Thus, dissipation is 
beyond the scope of the classical variational method [l| . 
To incorporate dissipative effects, we may introduce, for 
example, the Rayleigh dissipation function or some time- 
dependent external factor in actions, but it is not easy to 
specify uniquely them 0. 

Three decades ago, a new variational approach called 
the stochastic variational method (SVM) was introduced 
0, 13] • There, dynamical variables are extended to 
stochastic ones. In this way, dissipation is naturally in- 
duced by noise, without requiring additional functions. 
For example, it was shown that the Navier-Stokes (NS) 
equation for viscous fluids can be derived for an incom- 
pressible case d, @. SVM is also used to derive the 
Schrodinger equation from an action of classical dynam- 
ics. In the latter case, SVM has been discussed in rela- 
tion to Schrodinger and Nelson's formulation of quantum 
mechanics Q • If the general criteria of the applicability 
of SVM is established, it will serve as a very powerful 
tool in practice. This is because, as we will see later, 
once we know the form of Lagrangian of the reversible 
dynamics, the effect of dissipation to this dynamics is 
introduced in a systematic manner. In general, the con- 
struction of Lagrangians is relatively easier than that of 
evolution equations, because the form of the Lagrangian 
is constrained from the symmetry principle. 

In spite of these interesting concepts and features of 
the method, studies on its potential applicability to real- 
istic problems have not been sufficiently developed. For 
example, even for NS equation, only the case of an incom- 



pressible fluid has been studied. To claim the reliability 
of SVM, it is necessary to show explicitly that this can 
reproduce established dissipative equations. 

In this letter, we show that the full NS equation, at 
least, can be derived in the framework of SVM. Our ac- 
tion is the non-relativistic version of that used for the 
variational formulation of relativistic hydrodynamics 0] , 
but differs from those used in previous works of the clas- 
sical variational formulation of hydrodynamics [l|, Q . 

II. ACTION FOR IDEAL FLUID 

We first derive the Euler equation using the classical 
variational method. The same action used here is em- 
ployed for the derivation of the NS equation in SVM. 

Let us consider a fluid and divide it into small mass ele- 
ments specified by their Lagrangian coordinate R's. The 
space position of the fluid clement (the Euler coordinate) 
is denoted by 



r(R,i). 



(1) 



Let V = v(R, t) be the velocity field of the fluid element. 
When the time evolution is smooth, the time derivative 
of the Euler coordinate associated with the fluid element 
R is given by this velocity fleld (see later discussion). 



ar(R,t) 



dt 



r{-R,t). 
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Let U be the specific internal energy per unit mass of the 
fluid element. In general, in local thermal equilibrium, 
U depends on (R, t) only through the specific entropy s 
and the mass density p. With the energy density e, we 
can write U = ej p. The Lagrangian is then given by. 



i: = | rf^Rpo Qv2(R,t)-£^ 



(3) 



where the first and second terms represent, respectively, 
the kinetic energy and the "potential energy" associated 
with the fluid element. Here, po is the mass density mea- 
sured in the Lagrangian coordinate system and, by def- 
inition, does not depend on time. Note that the same 
Lagrangian can be expressed with the integral over the 
Eulcr coordinates as 



(4) 
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because, from Eq.(lT|), the two mass densities p and po 
are related with the Jacobian J = det|9r/9R| through 
the coordinate transformation as p = po/>^ 0- 
Then the action is expressed as 



dt 



(5) 



In the usual derivation of the Euler equation, the spe- 
cific entropy s should be kept constant for the variational 
procedure, representing that the fluid is ideal. Then we 
consider only the variation of r, leading to 



'at 



R 



1, 

P 



-VrP = 0, 



(6) 



where P is pressure defined by the thermodynamic rela- 
tion under the assumption of local thermal equilibrium. 



P 



d{l/p) 



(7) 



and the notation Vr represents the gradient with respect 
to the Euler coordinates r. In the following, the symbol 
V is used only for the gradient with respect to r, and 
hence the index is omitted. Equation ^ is the Euler 
equation, since 

9,v(R,i)|j,= [at+v(r,t).V]v(r,t). (8) 

In the above derivation, we used the following relations. 



d 



0, 



d 



where A^^ = dJ / d{dr'' / dlV). 



(9a) 
(9b) 
(9c) 



III. STOCHASTIC VARIATIONAL METHOD 

Following the spirit of SVM, we start from the same 
Lagrangian ([3]) as the ideal case. In SVM 0,0], however, 
we allow random fluctuations of r in its time evolution 
due to noise, so that the time derivative of r as Eq. ^ is 
not well-defined. The physical reason for the appearance 
of noise is the microscopic degrees of freedom which are 
coarse-grained in macroscopic scales. They act as the 
origin of fluctuations of the movement of fluid elements, 
leading to dissipation. Thus, the time evolution of r is 
described by the following stochastic differential equation 
(SDE), 



dr(t) = u rfi + B*dW(t). (dt > 0), 



(10) 



where, in the right hand side, the first term u is an un- 
known function of (r, i) which is to be determined by 
the variational procedure, and the last term is the noise 
term given by a Wiener process W(t) = {Wx,Wy,Wz) 
satisfying 



E[W,{s)-W,iu)] 
E [{W,{s) - W,{u)){Wk{s) ~ Wk{u))] 
E[{W,{s)~W,{u))Wk{t)] 



0, 
0, 



(11a) 
|s- 7.1, (lib) 
(11c) 



where k = x,y, z, — oo < t < s,u < oo and E [ ] de- 
notes the expectation value of stochastic processes. The 
symbol * is used to denote the Ito definition of a prod- 
uct for stochastic variables [l^l- The coefficient B is a 
vector composed of second rank tensors. For the sake of 
simplicity, we consider it as constant. We call Eq. ([TU]) 
the forward SDE since it is defined only for dt > 0. 

In SVM, the forward SDE is not sufficient to com- 
plete the formulation of stochastic variations. Because 
stochastic trajectories are not differentiable, the defini- 
tion of velocity is not unique. There are two possible 
definitions for the velocity of a fiuid element R at t, 

r{R,t + dt)~r{R,t) 
•vp — J> hm ; , 12a) 

dt^o+ dt 

V. ^ lim riR,t + dt)-riR,t)^ 

dt^o- dt 

When r is continuous and smooth, the two definitions 
should coincide, vj? = v^, as is the case of Eq. 
However, stochastic r is not smooth and we should dis- 
tinguish the two evolutions defined by v^^ and vg. 

The backward SDE which describes the time reversed 
process of Eq. ([TU]) is given by, 

dr(t) ^udt + Bi.dW{t), {dt<0) (13) 
where the new drift term u should be related to u, 



u 



u 



j 



2i/'^V^lnp. 



(14) 



This relation is obtained from the consistency condition 
of the two Fokker-Plank equations obtained from the two 
SDEs (Uni) and ((131) d, g. Here we introduced v'^ = 
[BB^]*^/2. The noise term W is again the Wiener pro- 
cess given by Eq. (1111) . There is no correlation between 
W and W. 

It should be emphasized that the velocities u and u 
are not parallel to the current of the mass density. From 
the Fokker-Plank equation obtained from Eq. (|10p . the 
mass density equation is uniquely given by 

dtP=-Y^ ^^{pu' - ^"'^oP) = -V • (pv„). (15) 



Here is parallel to the mass current, and defined by, 
v:„ = (u^ + u^)/2 = u' - ^ i/^d, In p. (16) 
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We call Vjn and u the mass velocity and the diffusion 
velocity, respectively. For an incompressible fluid, all of 
three velocities, u, u and Vm coincide due to Vp ~ 0, 
showing the intrinsic difference between the compressible 
and incompressible cases. 

Following Ref. Q , the action expressed with stochastic 
variables is then obtained from Eq. ([3]) by replacing v 
with the mean forward derivative Dr, 



dt 



lpr(R,t)).(i?r(R,i))-- 
2 p 



Here the mean forward derivative Dr is defined by 

'r{t + h)~r{t) 



Dr(t) = lim E 



(17) 
(18) 



Here, E[F{t')\Vt] denotes the conditional average of the 
time sequence {F {t') ,ta < t' < tf,}, taking the expecta- 
tion values for t' > t, fixing F {t') of t/ < t Q. Thus 
DF{t) is a stochastic variable. Note that the product of 
the mean forward derivatives is independent of the choice 
of the discretization scheme such as the Ito, Stratonovich- 
Fisk and Hanggi-Klimontovich schemes by the definition 
of the mean forward derivative [loj . 

The mean backward derivative is, similarly, defined as 



Dr(t) = lim E 



r{t) - r{t - h) 



h 



Ft 



(19) 



where, E[F{t')\Ft] is the conditional average for the past 
sequence, fixing the future values. 

There exist several ways to express the classical kinetic 
term in terms of Dr and Dr. For example, Dr-Dr is used 
in Ref. 0] . The final results for the equation of motion 
depend on this choice. To obtain the NS equation, Eq. 
([T7)) is employed. 

Similarly to the Euler equation, we consider variations 
for r as follows. 



rx{R,t)=r{R,t) + Xi{r,t), 



(20) 



where f{r,t) is an arbitrary function with the boundary 
condition, i{r,t)t=t^ = i{r,t)t=tt = and A is a smaU 
parameter. Then, for example, the variation of the mass 
density is calculated as 



6xp^ p\~p^ -X-jYl 



dJ 



didryOKJ) dW 



MX'). (21) 



Keeping the terms up to first order in A, d\I is given 



.6 

Sxl = - dtd^RY^ 

J a 



Po 



xE 



AP (DDr\R,t) + -4lP) + -SxS 
\ p or* / m 



, (22) 



where T = dmU /ds and m is the mass of molecules. In 
this derivation, we used the stochastic partial integra- 
tion formula, see Appendix. From Eq. (fT8|) and the Ito 
formula [lo| . we obtain 



DDr' = Du' =dtu' + Y^ dju' - 



v^'^djdku'. (23) 



Here the above expression corresponds for the general 
noise tensor B. For the following discussion to derive the 
NS equation, we consider the case of B = \/2i/I, where I 
is a unit matrix. 



IV. ENTROPY VARIATION 

Differently from the ideal case, s should be treated as 
a functional of r and we need to specify Sxs in terms of 
Eq. (j20p . For this purpose, we employ a following sim- 
ple model. Suppose that the fluid-dynamical time scale 
Thyd = p/p much larger than that of microscopic de- 
grees of freedom, Tmic, inside a fluid element. Then local 
thermal equilibrium should be achieved, recovering the 
ideal fluid where s is constant {6\s = 0). Therefore, we 
expect that the entropy variation is expressed in powers 
of Tmic/Thyd in such a way that the ideal fluid case is 
recovered in the vanishing limit of T,nic- We then write 

S\s = Sx{aiTjrdcp/p + a2{T„iicpl pY H )i where a/s are 

expansion coefficients. The lowest order truncation gives 



5\s = 5x {g{p)p) , 



(24) 



where g(p) is an arbitrary function of p. For the stochas- 
tic variation, p is interpreted as {D + D)p/2. Note that 
Sxs is the virtual change of s associated with the varia- 
tions and does not necessarily satisfy the thermodynamic 
principles such as Sxs > 0. 



V. NAVIER-STOKES EQUATION 

By substituting Eq. dM]) into Eq. (g^l), the condition 
of 5x1 = for the arbitrary function f (r, t) leads to 

p{dt + • V)u' + a,(P - /iV • v™) - ^ d,irid,u') = 0. 

(25) 

Here we used that p — — p'^ g{p) / ni{dT / d p) g and r] = up. 
The contribution from dxs effectively changes pressure 
by /iV • v,„. The coefficient p is known as the second 
coefficient of viscosity. 

As was pointed out, the fiuid velocity of the NS equa- 
tion is not u but v^. Eliminating u using Eq. (|16p . Eq. 
([25]) is finally re-expressed as 

Pidt + v™ • V)<„ + ^ d,[iP - CV • v™)<5^^- - ,?e^] 



J2d,(r^d,^jdanp))=0, 



(26) 
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where 



9j< 



(27) 



Wc thus identify rj as the shear viscosity and C = fi+2ij/3 
as the bulk viscosity. The last term is not only of second 
order for the magnitude of fluctuations v = rj/p, but also 
of third order for the spatial derivative V. In accordance 
with the approximation used in the NS equation, this 
term should be discarded as a higher order correction. 
In this sense, Eq. (pS)) is completely equivalent to the 
compressible NS equation. 



VI. CONCLUDING REMARKS 

In this letter, we showed for the first time that the 
stochastic variational method can be used to derive the 
full Navier-Stokes equation. The basic ingredients of this 
approach are the classical action of the ideal fluid and the 
stochastic motions of fluid elements which are induced by 
the white noise. 

The present derivation shows clearly that there are two 
different physical origins for the shear and bulk viscosi- 
ties. The shear stress tensor and a part of the bulk vis- 
cous pressure are obtained through the stochastic mo- 
tions of fluid elements, while the entropy variation in the 
potential term affects only the bulk viscous pressure (the 
second coefficient of viscosity). It is easy to see that the 
entropy variation due to the shear stress tensor e™ gives 
only higher order corrections. Thus, the shear viscosity 
in the NS equation is not altered by such a modification. 

So far, we have emphasized the mathematical aspects 
of SVM and did not discuss its physical background. To 
see it, note that the noise introduced in Eq. ((TU)) is re- 
lated directly to the transport coefficients. For example, 
we can show that rj satisfies the Einstein relation. 



P 



dtE[69{t) ■ 5v(0)], 



(28) 



where Sv ~ dr/dt — u =^/2i'dW /dt. This is nothing but 
the realization of the fiuctuation-dissipation theorem and 
it appears as a natural consequence of SVM. We thus 
conclude that SVM possesses not only the well-defined 
mathematical structure but also a reasonable mechanism 
of dissipation. 

Above results show that the SVM approach is consid- 
ered as a promising framework, and can be extended to 
more general dissipative phenomena which are not in the 
scope of the NS equation. In fac t, g eneralization of a 
diffusion equation is done in SVM . Another possible 
example can be found in the soft matter physics. One 



formulation of such dissipative equations is based on On- 
sager's variational method [l^l so that the comparison 
of SVM to Onsager's method will clarify the physics of 
variational approaches for the dissipative phenomena. 

As shown, the present result of SVM specifies the form 
of the higher order correction to the NS equation. If SVM 
is a reliable approach, this higher order term neglected 
in the NS equation should be considered seriously. The 
structure of Eq. ipS)) reminds that of the generalized hy- 
drodynamics proposed by Brenner [isl . but in the 
SVM approach, the difference of the two velocity fields 
appears as the higher order correction to the NS equa- 
tion. One example of the importance of higher oder terms 
obtained by SVM plays in fact a crucial role in diffusion 
processes. The application of SVM to a diffusion pro- 
cess leads to a generalized form of the diffusion equation 
[isj which contains memory effects. The existence of the 
higher order term guarantees Pick's law . 

As mentioned above, in spite of very attractive fea- 
tures, the general availability of SVM is not yet estab- 
lished. This question should be investigated by the ap- 
plications of SVM to field theoretical systems, relativistic 
systems and complex fiuids. 

This work was financially supported by CNPq, 
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Appendix A: stochastic partial integration formula 



The time variable is discretized as 

tj =a + j , j = 0,1,2, - •• , 

n 

Then we can show 

dtE[{DX{t)} ■ Y(t) + X(t) • DY(t)] 

n-l 



(Al) 



= lim E 

n 

+ \im E 



(X,+i-X,)^^^±i^ 



b — a 



n-l 



Xj + Xj-i 



(Y, - Y,_i) 



n 

b — a 
n 



= lim ^ i?[X,+i Y, - X, Y,_i] 



= E[X{b)Y{b) -X{a)Y{a)]. 



(A2) 



This is called the stochastic partial integration formula 

M. 
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